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ABSTRACT 


The Cosmic Evolution Early Release Science (CEERS) program observed the Extended Groth Strip 
with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In 
this paper, we discuss the four MIRI pointings that observed with longer wavelength filters, including 
F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with 
the Spitzer/MIPS 24 um population in the EGS field. We find that MIRI can observe an order of 
magnitude deeper than MIPS in significantly shorter integration times, attributable to JWSTs much 
larger aperture and MIRI’s improved sensitivity. MIRI is exceptionally good at finding faint (Lig « 
101? Le) galaxies at z ~ 1— 2. We find that a significant portion of MIRI galaxies are “mid-IR 


weak"-they have strong near-IR emission and relatively weaker mid-IR emission, and most of the 
star formation is unobscured. We present new IR templates that capture how the mid-IR to near-IR 
emission changes with increasing infrared luminosity. We present two color-color diagrams to separate 
mid-IR weak galaxies and active galactic nuclei (AGN) from dusty star-forming galaxies and find 
that these color diagrams are most effective when used in conjunction with each other. We present 
the first number counts of 10 um sources and find that there are S 10 IR AGN per MIRI pointing, 
possibly due to the difficulty of distinguishing AGN from intrinsically mid-IR weak galaxies (due to 
low metallicities or low dust content). We conclude that MIRI is most effective at observing moderate 
luminosity (Lig = 10? — 109 Le) galaxies at z = 1 — 2, and that photometry alone is not effective at 


identifying AGN within this faint population. 


1. INTRODUCTION 


In the past two decades, the Multi-Band Imaging Pho- 
tometer (MIPS) on the Spitzer Space Telescope revealed 
a new understanding of the infrared Universe through 
observations of a previously undetected population of 
infrared-luminous galaxies at z > 0.5 (Chary et al. 
2004; Papovich et al. 2004). MIPS observations demon- 
strated that the bulk of star formation in the Universe 
occurred when the Universe was a mere 3 — 6 billion 
years old (z ~ 1 — 2), the so-called “cosmic noon” epoch 
(Pérez-González et al. 2005; Caputi et al. 2007). Fur- 
thermore, contrary to the local Universe, this star for- 
mation occurred predominantly in dust rich galaxies- 
namely, (Ultra) Luminous Infrared Galaxies ((U)LIRGs, 
Lig > 10!! Lo, where Lig is the integrated luminosity 
from 8 — 1000 um; Magnelli et al. 2011). Although such 
galaxies exist today, they contribute relatively little to 
the local star formation rate (Rodighiero et al. 2010). 

While (U)LIRGs at cosmic noon are forming stars 
at prodigious rates (star formation rates, SFRs, of 
100 — 1000 M /yr), they predominately lie on the main 
sequence of star formation, in stark contrast to local 
(U)LIRGs (Elbaz et al. 2011; Kartaltepe et al. 2012). 
The main sequence is an empirical tight correlation 
between the SFR and stellar masses (M.) of galax- 
ies (Noeske et al. 2007; Elbaz et al. 2007). Galaxies 
on the main sequence are presumed to be undergoing 
long-lived secular evolution, while galaxies above the 
main sequence are undergoing short bursts of star for- 
mation, likely triggered by a major merger event. The 
Infrared Spectrograph (IRS) instrument on Spitzer was 
used to probe the mid-IR emission of these cosmic noon 
(U)LIRGs in detail and found that their mid-IR dust 
emission is remarkably similar to isolated LIRGs in the 
local Universe, indicating a metal-enriched environment 
and the distinct lack of an obscured nuclear starburst 
(Sajina et al. 2012; Kirkpatrick et al. 2012, 2015; Pope 
et al. 2013; McKinney et al. 2020). Simply put, dusty 
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star-forming galaxies in the distant Universe resemble 
scaled up versions of non-merging star-forming galaxies 
in the local Universe. 

All massive galaxies are presumed to host a supermas- 
sive black hole at their centers. The peak of the super- 
massive black hole accretion rate density also occurs at 
cosmic noon (Delvecchio et al. 2014; Peca et al. 2023), 
making this epoch of fundamental importance to under- 
standing the concurrent buildup of stellar and black hole 
mass. Furthermore, there is evidence from both deep 
X-ray and IR observations that the bulk of black hole 
growth at cosmic noon is obscured (Del Moro et al. 2016; 
Ananna et al. 2019; Yang et al. 2023a). Mid-IR emis- 
sion can reliably be used to identify both obscured and 
unobscured active galactic nuclei (AGN; Padovani et al. 
2017; Yang et al. 2023a). In previous mid-IR surveys 
with Spitzer/MIPS or the Wide-Field Survey Infrared 
Explorer (WISE), AGN make up at least 15% of extra- 
galactic sources with S24 > 100uJy (Jarrett et al. 2011; 
Kirkpatrick et al. 2017a, hereafter K17). 


In star-forming galaxies, the mid-IR emission is dom- 
inated by polycyclic aromatic hydrocarbon (PAH) com- 
plexes at 6.2, 7.7, 11.2, and 12.7 um. AGN typically have 
a dusty toroidal structure surrounding the accretion disk 
that radiates as a power-law in the mid-IR. These dif- 
ferences make the mid-IR advantageous for identifying 
populations of star-forming galaxies and AGN, and sev- 
eral color-color methods exist for doing so (Lacy et al. 
2004; Stern et al. 2005; Donley et al. 2012; Kirkpatrick 
et al. 2013; Messias et al. 2012). The Mid-Infrared In- 
strument (MIRI) (Rieke et al. 2015) on the James Webb 
Space Telescope (JWST ; Gardner et al. 2006) has several 
improvements over MIPS: the greatly improved spatial 
resolution allows for dramatically improved source de- 
blending, and the MIRI filter set better spans the mid-IR 
emission. These improvements led to studies predicting 
that MIRI would be efficient and reliable at finding AGN 
(Messias et al. 2012; Langeroodi & Hjorth 2023) at high 
redshift (z > 5; Volonteri et al. 2017) or at lower lumi- 
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Table 1. CEERS Red MIRI Pointings 


Pointing i RA Dec Filter Exp. Time 50 depth b 

(32000) (32000) (sec) (uJy) 

1 14:20:38.88 +53:03:04.6 770W 1648 0.21 
1000W 1673 0.42 

1280W 1673 0.78 

1500W 1673 1.27 

1800W 1698 2.51 

2100W 4812 4.79 

2 14:20:17.42 +52:59:16.2 770W 1648 0.21 
1000W 1673 0.43 

1280W 1673 0.82 

1500W 1673 1.15 

1800W 1698 2.73 

2100W 4812 3.87 

5 14:19:03.74  4-52:49:08.5 1000W 1243 0.50 
1280W 932 1.27 

1500W 932 2.31 

1800W 1243 3.98 

8 14:19:21.78 +52:48:55.5 1000W 1243 0.47 
1280W 932 1.32 

1500W 932 2.00 

1800W 1243 3.98 


“For layout of the pointings, see https://ceers.github.io/obs.html 


b Flux limit calculations are derived from the AB magnitude limits presented in Yang et al. (2023b) 


nosities and high obscuration levels (Kirkpatrick et al. 
2017a; Yang et al. 2021). 

In this paper, we explore the nature of the sources in 
one of the first MIRI surveys, which was observed as part 
of the Cosmic Evolution Early Release Science (CEERS; 
PIS. Finkelstein) program. In Section 2, we discuss the 
observations. In Section 3, we compare the MIRI sources 
to the MIPS sources in the Extended Groth Strip (EGS) 
field, and we discuss the presence of AGN. We discuss 
the implications of the MIRI observations, including the 
surprising lack of AGN, in Section 4. Finally, we present 
our conclusions in Section 5. 

In this work, we assume a standard flat cosmology 
with Ho = 70 Mpc/km/s, Qm = 0.3, and A = 0.7. 


2. DATA 
2.1. MIRI Observations 


The CEERS MIRI observations, data reduction, and 
photometric catalogs are fully described in Yang et al. 
(2023a); here, we summarize relevant details. The EGS 


field was observed with eight MIRI pointings. Four 
“blue” pointings were only observed in the F560W and 
F770W filters. In this paper, we discuss the four “red” 
pointings, observed in the longer wavelength filters (Ta- 
ble 1). Pointings 1 (Figure 1) and 2 were observed 
in July 2022, while pointings 5 and 8 were observed 
in December 2022. Pointings 1 and 2 do not over- 
lap with CEERS NIRCam observations. The F770W 
and F2100W filter are only observed in pointings 1 and 
2. Table 1 lists the 5c flux limits for each filter in 
each pointing (Yang et al. 2023b). The MIRI FOV 
is 112.6” x 73.5”, bringing the total observed area to 
9.2arcmin? for the 4 MIRI pointings discussed in this 
work. 

The MIRI photometry was extracted with T-PHOT 
(Merlin et al. 2015) using Hubble Space Telescope 
(HST)/F160W observations as a prior. HST positions 
were used rather than NIRCam due to the lack of uni- 
form NIRCam coverage in the MIRI pointings. For de- 
tails of the MIRI observations, data reduction, photom- 
etry extraction, and quality assessment see Yang et al. 


53:02:00.0 


30.0 


34.0 


Figure 1. We show MIRI pointing 1 (right panel) alongside the Spitzer/IRAC (middle) and MIPS (left) observations of the 
same region. The apertures show the location of detected sources in each image (MIRI region only). For the MIPS (IRAC) 
image, the apertures are 6” (2”), corresponding to the instrument beam size. In the IRAC image, blue corresponds to channel 
1 (3.6 um), green corresponds to channel 2 (4.5 jum), and red corresponds to channel 3 (5.8 um). In the MIRI image, the 770W 


filter is blue, F1000W is green, and F1280W is red. 


(2021); Papovich et al. (2022) and the CEERS MIRI 
observation paper (Yang et al. 2023b). 

The use of near-IR HST priors resulted in some spu- 
rious detections where a galaxy is visible in the HST 
image but not in all of the MIRI filters. To be included 
in the final sample, we required sources to 1) have a SNR 
> 3 in at least two MIRI bands and 2) be well-fit by a 
mid-IR template (see Sections 2.4 and 3.2 for a full de- 
scription of the template fitting procedure). There are 
911 sources with SNR > 3 in at least one MIRI band. 
Of these, there are 575 (466) that have SNR » 3 in at 
least two (three) bands. From the parent sample of 575 
sources, 482 were well-fit by mid-IR templates and are 
included in the final sample. The remaining 93 sources 
have faint fluxes (< lyJy) in the mid-IR (A > 5m) 
and no unambiguous PAH features. 


2.2. Redshifts 


We cross-match our MIRI sources with the CAN- 
DELS (Grogin et al. 2011; Koekemoer et al. 2011) 
EGS catalog of Finkelstein et al. (2022), which in- 
cludes HST/ACS+WFC3 and Spitzer/TRAC photom- 
etry and photometric redshifts derived with the EAZY 
code (Brammer et al. 2008), fitting a suite of tem- 
plates to the HST+IRAC photometry. The code pro- 
duces a redshift posterior probability distribution which 
includes the 68% confidence intervals. The final red- 
shift is the peak of the probability distribution. As dis- 
cussed below, for some MIRI sources we use photomet- 
ric redshifts from the CANDELS catalog of Stefanon 
et al. (2017), measured by fitting HST and IRAC data 
with various different codes. The Zphot is the median 
redshift determined through five separate codes that fit 


templates to the UV/optical/near-IR data. A mere 20 
sources have spectroscopic redshifts, and we use these 
where available. 

Optical photometric redshifts can carry a degree of 
uncertainty if templates are not representative of the 
full galaxy population or if the photomery is noisy or 
sparse (Newman & Gruen 2022). MIRI covers the PAH 
features for the majority of our galaxies, which may pro- 
vide additional constraints on the redshifts due to their 
prominence (Pope et al. 2008; Kirkpatrick et al. 2012). 
MIRI is the first instrument with closely spaced filters 
allowing for the identification of PAH features with pho- 
tometry alone. 

To test how well fitting templates to mid-IR data can 
reproduce optical redshifts, we fit all galaxies with three 
templates described in more detail in Section 2.4. Two 
templates are a star-forming template ( MIR0.3') and an 
AGN template (‘MIRO.7’) from the Kirkpatrick et al. 
(2015, hereafter, K15) library. The third template is 
the Lig = 10!? Lo template from the Chary & Elbaz 
(2001, hereafter, CEO01) library. We include the IRAC 
4.5 um and 5.8 um photometric points, along with all 
available MIRI data, in the fitting as these are useful 
in identifying the stellar bump at A ~ 1.6 jum and the 
stellar minimum at Á ~ 5ym. For all but 32 sources, 
the MIRI redshift agrees with the optical redshift from 
Finkelstein et al. (2022) within the confidence interval 
on the optical redshift. We therefore adopt the Finkel- 
stein et al. (2022) redshifts for these galaxies. Of the 32 
sources where the two redshifts do not agree, the MIRI 
redshifts of 24 sources agree with the optical redshifts 
from Stefanon et al. (2017) within the confidence inter- 
val on the optical redshift. We adopt the Stefanon et al. 


CEERS KEv PAPER VII 5 


(2017) redshifts for these sources. For the remaining 8 
sources, we adopt the MIRI redshifts. 


2.3. MIPS 24 um Observations 


The EGS field was observed with Spitzer/MIPS 
as part of the Spitzer/FIDEL Legacy project! (PI 
M. Dickinson; Magnelli et al. 2009). The typical 
MIPS exposure time per position within the CAN- 
DELS/CEERS areas is 14,000s, and the MIPS data 
analyzed here cover the CANDELS WFC3 field of 
„200 arcmin?, nearly 20x larger than the CEERS MIRI 
area. The MIPS beam size is 6", causing blending of 
closely spaced sources within each beam, which have 
been carefully deblended based on HST and IR AC prior 
positions and fluxes. Photometry has been calculated 
with a PSF-fitting method to extract MIPS fluxes at 
the positions of the HST/F160W priors from the cata- 
log of Stefanon et al. (2017). The methodology is similar 
to what is described in Jin et al. (2018) and Liu et al. 
(2018), and full details will be provided in a forthcom- 
ing paper (Le Bail et al.in preparation). Flux biases 
and uncertainties were measured using extensive simula- 
tions. There are „4000 MIPS sources with S24 > 21yJy, 
which is the average 3c flux limit in the Le Bail catalog. 

Figure 1 illustrates the image quality of MIRI com- 
pared with IRAC and MIPS. With MIPS, the larger an- 
gular resolution leads to multiple galaxies falling within 
one beam. MIRI observations alleviate this challenge, 
enabling the detection of fainter, non-isolated sources. 


2.4. Infrared Templates 


To measure redshifts, estimate Lyr, and search for 
AGN emission, we fit a suite of templates to the MIRI 
population (Section 3.2). These templates are the K15 
MIR Library and the CEO1 Library. We opt for this 
method as it is less computationally expensive than 
codes such as CIGALE. The K15 and CE01 templates are 
empirically derived from infrared observations of galax- 
ies, while CIGALE and other spectral energy distribution 
(SED) fitters produce model spectra using stellar popu- 
lation synthesis, dust emission models, radiative transfer 
calculations, or energy balance requirements. 

The K15 templates were empirically created from 343 
24um-selected galaxies spanning the redshift range z = 
0.5—3.0. Every galaxy had a Spitzer/IRS spectrum (low 
resolution spectroscopy) that was modeled by combining 
a PAH template with power-law emission attributed to a 
torus. This spectral modeling determined f(AGN)wrn, 
which is the fraction of mid-IR emission attributable 
to the power-law component. f(AGN) ir spans the 
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range 0.0-1.0 (fully star-forming to fully AGN). Addi- 
tionally, all galaxies had JHK, Spitzer/IRAC, and Her- 
schel/ PACS and SPIRE photometry covering the full IR 
wavelength range. 

Galaxies were sorted by f( AGN)wrn in increments of 
0.1 (that is, the pure star-forming bin would have galax- 
ies with OX f(AGN)yrm «0.1, while the pure AGN bin 
would have 0.9 f( AGN)yin €1.0; every other bin con- 
tains galaxies with f(AGN) ir indicating varying mixes 
of star-formation and AGN emission) to create empirical 
templates. From A = 0.9 — 20 um, the normalized pho- 
tometry and spectroscopy were averaged together as a 
function of wavelength using a bootstrapping technique. 
Beyond 20 um, the normalized photometry was fit with 
a two-temperature modified blackbody. 

Crucially, all galaxies in the sample used to create 
the K15 template library have Lig > 10! Lọ. They 
are representative of the (U)LIRG population at cosmic 
noon. 

As discussed below, lower luminosity templates are 
also needed to represent the MIRI population. To this 
end, we included the CE01 templates with Lig = [2 x 
108,5 x 105,1 x 109,5 x 109,1 x 101°, 5 x 101°] Le in our 
fitting procedure. 

CEO01 used models tailored to match the ultraviolet- 
submillimeter SED of four galaxy prototypes: Arp 220, 
NGC 6090, M82, and M51. These four representative 
models were then divided into mid-infrared (4-20 um) 
and far-infrared (20-1000 um) components and fit to the 
ISOCAM and IRAS observations of local galaxies and 
(U)LIRGs in the IRAS Bright Galaxy Survey. Galax- 
ies were sorted into luminosity bins, and the best-fitting 
mid-infrared and far-infrared models were averaged to- 
gether to yield the final representative template SED for 
each luminosity bin. 


3. RESULTS 
3.1. Comparison with Spitzer 24 um 


Prior to the launch of JWST, Spitzer was the pre- 
mier mid-IR observatory. To better understand the new 
mid-IR Universe revealed through MIRI observations, 
we compare the properties of the CEERS JWST/MIRI 
galaxies (traced through the 18 um population) to the 
EGS Spitzer/MIPS population. 

Figure 2 compares the luminosity and redshift distri- 
bution of the full EGS MIPS population and the 18 jum 
population. We choose to show the 18 um population 
as 18 um is close in wavelength to 24 um, F1800W was 
observed in all 4 pointings, and the F1800W sensitivity 
is better than at F2100W. There are 433 sources with 
SNR> 3 in the F1800W filter. 
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Figure 2. Redshift vs. vL, for the full MIPS population (v = 24 um; black crosses) and the MIRI popoulation (v = 18 um; 
orange squares). The surprising lack of high-z MIRI sources is attributable to the lack of bright dust features at Arest < 5um 
and the much smaller survey area. MIRI is capable of detecting fainter sources than MIPS at all redshifts in a fraction of the 
time. We show the final sample (SNR 3 in at least 2 bands and well-fit by a mid-IR template; 385 18 jum galaxies meet this 
criteria) in orange. We also show the sources which have SNR>3 in only the F18000W filter (14 sources; blue circles with 
crosses), and the sources which are not well-fit with a mid-IR template (34 sources; blue triangles). Our selection criteria is not 
biasing the final sample towards brighter sources. The histograms in the side panels have been arbitrarily normalized due to 


the stark difference in numbers of sources. 


As seen in the top panel, MIPS sources have a bi- 
modal redshift distribution with peaks at z ^ 1 and 
z ~ 2 due to PAH emission. At z ~ 1, the 11.2 4m and 
12.7 um PAH complexes fall within the 24 um bandpass, 
and at z ~ 2, the 7.7 um PAH feature slides in. Addi- 
tionally, the 9.7 jum silicate absorption feature falls in 
the bandpass at z ~ 1.5, contributing to the dearth of 
detections. The 18 um distribution also has noticeable 
peaks where different PAH features slide into the filter. 
The peak around z ~ 0.5 is attributable to the 11.2 
and 12.7 um features. The majority of 18m sources 
lie around z ~ 1.3. This is due to the 7.7 um complex 
falling in the F1800W bandpass at z = 1.3. The 7.7 um 
complex is the broadest of the PAH features, spanning 
A = 7.2— 8.2 um. This rest wavelength range will cause 


a brightening at Aops = 18 jum from z = 1.2 — 1.5. The 
6.2 um PAH feature falls in the F1800W bandpass at 
z — 1.9, also visible in the redshift distribution. 

'The redshift distributions are most instructive when 
considered alongside the vL, distributions. Here, vL, 
refers to the restframe photometry in either the 24 um 
or F18000W filter. MIRI galaxies are an order of magni- 
tude less luminous in the infrared, at all redshifts, which 
can be attributed to the increased sensitivity. It is in- 
teresting, and perhaps unexpected, to note the lack of 
MIRI detections of bright sources at z > 3. At z > 3, 
F1800W probes the faintest dust emission (A S; 5um), 
contributing to the lack of high redshift sources. In 
terms of numbers of high—z sources, MIPS has the ad- 
vantage of a much larger FOV than MIRI (12x greater). 
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Uncovering a high redshift population of dusty galaxies 
will require a large area MIRI survey with the longest 
wavelength filters in order to probe the dust emission at 
higher redshifts. 

The 18 um population has a very similar luminosity 
and redshift distribution as the 21 wm population. There 
are 121 sources in pointings 1 and 2 detected with MIPS 
and 264 sources detected in MIRI F2100W. The 21 um 
MIRI distribution peaks at more than an order of magni- 
tude lower flux than the MIPS distribution, attributable 
to the improved sensitivities of JWST compared with 
Spitzer, largely due to the nearly 8x increase in diam- 
eter. F1800W reaches fainter flux limits than F2100W 
in similar exposure times, but F2100W can trace dust 
emission out to higher redshifts. 

Figure 2 highlights MIRI’s strengths: low luminosity 
galaxies at cosmic noon (z = 1 — 2). Here, MIRI probes 
a population inaccessible from Spitzer, and it does it in 
less than an hour. The lack of filters at Á > 21 ym mean 
that the brightest dust emission, the PAH features, are 
no longer detectable at z > 2. This makes MIRI ideally 
suited for detecting galaxies well below the knee of the 
luminosity function at cosmic noon. 


3.2. Determining What Powers the Mid-IR Emission 


Kirkpatrick et al. (2017a) predicted that MIRI ob- 
servations alone could be used to quantify AGN emis- 
sion. We test that prediction by estimating AGN frac- 
tion, f(AGN) ir, by initially fitting only the template 
library from K15 to the MIRI data. Template fitting is 
a complementary technique to full SED decomposition. 
The advantages of fitting templates are that it is less 
computationally expensive, there are fewer parameters 
resulting in fewer degeneracies, and it can be performed 
when only a few photometric observations are available. 

We fit the suite of templates to all galaxies that have 
two or more MIRI photometric points with SNR > 3.0. 
We only include the MIRI data in this initial fitting, 
since this is the range that the templates are most sensi- 
tive to. We only fit one parameter, which is the scaling of 
the template to the photometry. To assess the goodness- 
of-fit, we calculate a reduced y? parameter, and we ac- 
cept all fits with reduced x? = 0.5 — 2.0. We then vi- 
sually inspect the fits and overplot the Spitzer/IRAC, 
MIPS, and HST/WFC3 (F125W, F160W) photometry. 
98% galaxies with Lig > 10!! Lọ and 70% of galax- 
ies with 101? Lọ < Lir < 10! Lg are well-fit with 
the K15 templates. This is not surprising since the 
templates were created from IR-luminous galaxies with 
Lig ~ 1011 — 107? Lo. Compellingly, upon visual inspec- 
tion, the non-MIRI photometry follows the best-fitting 
template as well. 


However, the low luminosity galaxies probed by MIRI 
are very poorly described by the K15 template library 
below the mid-IR regime. Even when the reduced x? 
is acceptable, visual inspection shows that the observed 
near-IR emission is an order of magnitude brighter or 
more than the template emission. 

In order to better capture the emission of lower- 
luminosity galaxies, we expanded the template fitting 
to include templates from the CE01 Library. We fit 
the K15 and CE01 templates, this time including the 
Spitzer/IRAC 4.5 um data point in the fitting to help 
anchor the ratio of near-IR to mid-IR emission. We in- 
tegrate each best-fit template from 8 — 1000 um to esti- 
mate the galaxy's infrared luminosity. As stated in Sec- 
tion 2.1, we are able to achieve good template fits for 482 
galaxies, and these galaxies comprise the final sample. 
Interestingly, we find that the best-fitting CEO1 tem- 
plate Lyr does not necessarily correlate with the galaxy 
Lig derived from the renormalized SED templates. In 
~ 40% of cases, these numbers differed by more than an 
order of magnitude. 

In the initial fitting, where only the MIRI data points 
were fit with the K15 templates, a surprising 4396 of 
galaxies with Lig < 10! Lọ were classified as AGN 
(f(AGN)um 20.5). When the 4.5 um data point was 
included, along with the CE01 templates, this number 
dropped to 1096. Through our template fitting, we iden- 
tify 28 AGN (f(AGN)mrr>0.5) out of our final sample 
of 482 galaxies. For context, Yang et al. (2023a) fit UV- 
MIRI photometry with CIGALE for our same sample, and 
they identify 31 AGN. We compare our template fitting 
results with the CIGALE decomposition and find good 
overall agreement. For 6096 of sources, the AGN frac- 
tions agree to within 10%, and for 95% of sources, the 
fractions agree to within 5096. Broadly speaking, both 
manners of measuring f(AGN)mir will classify the same 
galaxies as AGN, although the exact fractions may vary. 

'The K15 templates were created from dust-rich galax- 
ies with strong PAH emission, and the ratio of mid-IR 
emission to near-IR emission reflects the fact that most 
of the galaxy's energy is being emitted in the infrared. 
In contrast, the lower-Iuminosity CEO01 templates repre- 
sents galaxies with less mid-IR emission compared with 
the strength of their near-IR emission. We label sources 
best fit with the CEO1 templates as “mid-IR weak" 
galaxies. Figure 3 shows how the AGN fraction and 
mid-IR weak fraction varies with luminosity and red- 
shift. As expected, the mid-IR weak fraction increases 
sharply below Lig = 10!? Lo. The fraction of mid-IR 
weak galaxies decreases with redshift, although this is 
likely driven by the flux limit of the survey. The fraction 
of galaxies that are AGN increases with redshift. The 
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Figure 3. The fraction of galaxies that are mid-IR weak (i.e, fit with CEO1 templates rather than K15; blue squares) or AGN 
(f(AGN)mir> 0.5; red triangles) as a function of Lig (Left) and redshift (Right). The small number statistics make it difficult 
to narrow down the true distribution of AGN with redshift, although it seems to be increasing. 


increase in AGN is similar to what was found in Yang 
et al. (2023a). The increase in the number of detected 
AGN is likely due to their enhanced emission around 
5um compared with star forming galaxies, making them 
easier to detect in the longer wavelength MIRI filters at 
higher redshift. 


3.3. Average Emission of Mid-IR Selected Galaxies 


MIRI has opened a deeper window on the infrared 
Universe, particularly at cosmic noon. We explore the 
ratio of near-IR to mid-IR emission in MIRI galaxies 
through the creation of average SEDs. We sort galaxies 
into luminosity bins of log Lig = 0.5dex, which allows 
for at least 35 galaxies in each bin. We then calculate 
the median [yp in each bin and normalize each galaxy 
to that value. The [yrs used to sort and normalize the 
galaxies are those derived by fitting the K15 and CE01 
templates. We discuss Lig further below. We remove all 
galaxies with f(AGN)ir> 0.3. We list the properties 
of each luminosity bin in Table 2. 


We include the HST/ACS (F606W, F814W), WFC3 
(F125W, F160W), Spitzer/IRAC, MIPS 244m, and 
MIRI photometry, redshifted to the rest-frame, in the 
creation of the SEDs. We use a bootstrap technique 
where we resample the sources in each bin with replace- 
ment 1000 times. In each iteration, we calculate the me- 
dian luminosity in wavelength bins of differential sizes. 
After 1000 iterations, we calculate the average luminos- 
ity in each bin (this is the average of the median lu- 
minosities calculated in each iteration) and the stan- 
dard deviation. We list these values in Table 3. Figure 


4 shows all of the galaxies in each luminosity bin and 
the resulting average SED. It is notable that the noisi- 
est photometry in each bin occurs around A ~ 5.0 um, 
where the emission switches from being dominated by 
the older stellar population to being dominated by dust. 
The highest luminosity bin shows the most uniformity 
among individual galaxies. 

Figure 5 shows the average SEDs in the left panel. 
The lower luminosity SEDs clearly have more of their 
energy being radiated in the near-IR, which could indi- 
cate that the bulk of their star formation is unobscured. 
We explore this further in Section 4. We test whether 
the shape of the SED changes from z = 0 to cosmic 
noon by comparing with the CE01 templates in the right 
panel of Figure 5. We plot the CE01 template with Lim 
closest to each average SED. Each CE01 template and 
average SED of a given luminosity have been arbitrarily 
normalized (the same normalization is applied to each 
set of a given luminosity) to allow for easier comparison. 
'The average SEDs are in general remarkably similar to 
the CEO1 SEDs, indicating that the relative amounts 
of near- and mid-IR emission does not evolve strongly 
with redshift for galaxies of the same infrared luminos- 
ity. However, it is notable that the shape of the near- 
IR emission does change quite strongly for most of the 
luminosity bins. The MIRI templates lack the strong 
emission seen at 1.6 um in the CE01 templates. This 
may reflect a difference in the ages of the stellar pop- 
ulations. The CEO01 templates were created from local 
galaxies which have had an additional = 6 — 9 billion 
years to evolve, which can make the stellar population 
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Table 2. Properties of Luminosity Bins and SEDs 


Bin Sources Median log Li ^ Median z Average SED log Lir? 
(log Lig) (La) (La) 
<9.0 50 8.52 [8.22,8.77] 0.58 [0.41,0.68 8.44 [0.55] 
9.0-9.5 54 9.26 [9.13,9.38] 0.90 [0.74,1.23 9.02 [0.06] 
9.5-10.0 79 9.74 [9.65,9.86] 1.33 [0.92,1.55 9.53 [0.07] 
10.0-10.5 77 10.24 [10.10,10.38] 1.32 [0.79,1.82 9.88 [0.06] 
10.5-11.0 43 10.72 [10.59,10.90] 1.42 [1.08,1.89 10.29 [0.04] 
11.0-11.5 39 11.14 [11.09,11.31] 1.57 [1.19,2.15 10.66 [0.04] 
> 11.5 35 11.75 [11.55,11.94] 1.53 [1.32,1.81 11.26 [0.04] 
@ Median log Lim and redshift [with lower, upper quartiles] for sources per bin. Lrr is 


extrapolated for individual sources using the best-fitting K15 or CE01 template. 


b Average Lin [with standard deviation] of log Lig values extrapolated by fitting the 
average SED per luminosity bin using several different template libraries. 


much older and the resulting emission much stronger in 
the near-IR. 

As a caveat—- Lig for individual sources was calculated 
by simply integrating the best-fit K15 or CE01 template 
from 8 — 1000um, and the Lyr for each bin listed in Ta- 
ble 2 is the median of these values. Very few galaxies 
in these CEERS MIRI fields have direct far-IR. detec- 
tions from Herschel, so Lig is an extrapolation from the 
mid-IR. If the ratio of mid-IR to far-IR emission is insen- 
sitive to redshift, then the Lig of fainter galaxies (those 
missed by Herschel) can be accurately estimated from 
local templates. 

We also calculate Lyp by fitting a suite of IR templates 
from CE01, Dale & Helou (2002), Draine & Li (2007), 
and Rieke et al. (2009) to each average SED (Table 3) 
at Arest > 4.7 um. Table 2 lists the average and stan- 
dard deviation of the log( Lig) values extrapolated using 
these four template libraries. The average log(Lrm) is 
about 0.5 dex smaller than the median log(Lig) values 
for individual sources per bin that were derived using the 
K15--CEO01 templates. Except for the lowest luminosity 
bin, the standard deviations are small (€ 0.1 dex), show- 
ing that existing SED libraries have broadly consistent 
mid-IR to total IR luminosity ratios, i.e., that the choice 
among current-generation SED templates does not have 
a large impact on total IR luminosities that we infer 
from MIRI data. 

Although the various template libraries largely agree 
with one another when used to extrapolate Lig from 
the mid-IR emission, it is important to bear in mind 
that these libraries have been calibrated on local sources 
or infrared luminous galaxies. The ratio of mid-IR to 
far-IR emission can change with location on the main 


sequence, redshift, and metallicity (Elbaz et al. 2011; 
Kirkpatrick et al. 2017b; Schreiber et al. 2018). 

Direct far-IR measurements for the low-luminosity, 
high redshift MIRI galaxies would be necessary to verify 
their total luminosities. While deep millimeter observa- 
tions with ALMA, LMT/TolTEC, or other facilities can 
provide some constraint at wavelengths beyond the peak 
of bolometric emission, ultimately a new, sensitive far- 
IR observatory will be required to robustly measure Lim 
for these galaxies. 


3.4. Color Selection 


Kirkpatrick et al. (2017a) used the K15 templates to 
predict where AGN would lie when combining MIRI 
photometric bands for color selection. MIRI color selec- 
tion of AGN is hampered by two facts: first, with just 
a few closely spaced photometric points, star-forming 
galaxies can mimic a power-law in the mid-IR. There- 
fore, MIRI color selection of AGN only has a unique so- 
lution when the redshift is known; second, if only mid-IR. 
photometry is used, mid-IR weak galaxies can mimic the 
weak PAH features in AGN. We find that when plotting 
the real CEERS observations in the MIRI colorspaces 
presented in Kirkpatrick et al. (2017a), the AGN regions 
are severely contaminated by the mid-IR weak galaxies. 

Instead, we present a new diagnostic that does not 
take into account redshift information, in order to 
illustrate the limited utility of using mid-IR colors 
alone to identify AGN. Figure 6 combines the colors 
Sr1800w / Sr1000w. With Sp1280w /.Sr1000w. We find that 
combining three filters, rather than four, does a better 
job at separating AGN, star-forming galaxies, and mid- 
IR weak galaxies. 
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Figure 4. To create average SEDs (shown in black), we separated galaxies by luminosity and normalized to the median 
luminosity in each bin. Individual galaxies are plotted by color. 


The dust-rich galaxies (those fitted with K15 tem- We define an irregular hexagon to separate the AGN 
plates) are plotted as the shaded circles, while the mid- from the non-AGN, the vertices of which are: 
IR weak galaxies (fit with CE01 templates) are plotted 


as the grey squares. AGN (solid circles) lie towards the 2 — log(Srisoow /$r1000w) 


center of the colorspace, while star forming galaxies lie y = log(Srizsow / Sr1000w) 
around the edges. [x, y] =[0.04, 0.15], [0.04, 0.44], 
[0.24, 0.79], (0.47, 0.79], 
[0.25, 0.44], [0.25, 0.15] (1) 


Within this region, there are 19 AGN and 48 contami- 
nants (there are 25 AGN total in this figure). Without 
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Figure 5. Left panel-Average SEDs from Figure 4. The PAH features clearly become stronger with increasing luminosity. 
At low-luminosities, the mid-IR emission is much weaker than the near-IR emission. Right panel-Average SEDs are plotted 
alongside the CE01 template of the same Lig. Arbitrary normalization is applied to templates of the same luminosities. The 
CE01 templates have more emission around A ~ 1.6 jum, which may be attributable to their older stellar populations. 


including redshift information, this reliability (2896) is 
the best that can be achieved with MIRI colors alone. 

Galaxies with PAH features move in the diagram as 
redshift increases. The star-forming galaxies move from 
the bottom of the plot at z « 0.6 to the top right at 
z ~ 0.6 — 1.5 (redshift track is illustrated with the pur- 
ple dashed line). At greater redshifts, they move to- 
wards the middle left. The mid-IR weak galaxies (red- 
shift tracks are the dot-dashed magenta and pink lines) 
follow the same general movement, but at lower color 
ratios, causing them to contaminate the region predom- 
inantly occupied by AGN. The AGN, due to their lack 
of PAH features, do not strongly change location due to 
redshift. 

As seen in the left panel of Figure 6, mid-IR weak 
galaxies contaminate the AGN region mostly at z ~ 1. 
The reason for this is illustrated in the right panel. 
We have plotted the K15 MIR1.0 template and our 
log Lig = 9.0 — 9.5 and log Lir = 10.0 — 10.5 SEDs from 
this work. We overplot the MIRI F1000W, F1280W, 
and F1800W filters, shifted to the rest wavelength they 
would have at z — 1. At this redshift, there simply 
is not a wide separation between the SEDs, so intrin- 
sic variation between sources will cause overlap between 
AGN and mid-IR weak colors. It is also obvious that 
adding a power-law criterion (i.e. F1000W < F1280W 
< F18000W) will not help the situation as that can be 
true for galaxies with and without PAH features in this 


regime. The part of the SEDs with the largest separa- 
tion is the near-IR. 

Figure 7 shows the filters where the three templates 
(plotted on the figure) have the strongest separation over 
the largest range of redshifts. The inclusion of a near-IR 
data point can help separate mid-IR weak galaxies from 
dust-rich galaxies. We combine IRAC and MIRI data 
to measure log Sr1800w / $3.6 m and log SP1000w / $3.6. 
We use IRAC rather than NIRCam due to the uneven 
coverage of the CEERS observations. 

The mid-IR weak galaxies (grey squares) lie in the 
lower left region of this colorspace, while the dust- 
rich galaxies lie in the upper left. The threshold 
Sp1000w/S3.6 > 0.35 separates dust-rich star-forming 
galaxies from the mid-IR weak galaxies. Given the dif- 
ference in luminosities between these classes of galaxies, 
it can also be seen as an infrared luminosity threshold. 
Dust-rich star-forming galaxies at z « 1 are the main 
contaminant in upper right, where the AGN predomi- 
nately lie. We recommend the following criteria for se- 
lecting AGN: 


x = log(Srisoow / $3.6) 

y = log(Smiooow / $3.6) 

y > 0.55 

y > 2. 78x — 0.98 

y € 1.92% + 0.45 (2) 
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Figure 6. Left panel MIRI color selection of AGN. Dust-rich galaxies (circles) are shaded according f(AGN)mir. Galaxies with 
f(AGN)mir> 0.5 are shown as the filled circles. The mid-IR, weak galaxies are the grey squares. The proposed AGN selection 
region (Equation 1) is shown as the solid line. We illustrate how the dust-rich star forming galaxies (purple dashed line, K15 
template) and mid-IR weak galaxies (this work, pink and light pink dot-dashed lines) move with redshift. AGN do not strongly 
evolve with redshift. Dust rich galaxies separate by f(AGN)mir, while the mid-IR weak galaxies contaminate the AGN regions. 
Right panel-mid-IR weak SEDs compared with a K15 AGN template. The IRAC Ch 1 (3.6 jum) bandpass is shown alongside 
the MIRI bandbpasses. All bandpasses have been redshifted to z = 1. At this redshift, the mid-IR weak and AGN templates 
have similar mid-IR colors, but the inclusion of a near-IR data point will help to separate them. 


Within this region, there are 15 AGN and 13 contami- 
nants. There are 22 AGN total shown in Figure 7. 

The MIRI-only selection (Figure 6) has a reliability 
of 28%, while the MIR+NIR selection (Figure 7) has 
a reliability of 5496. If we combine the mid-IR criteria 
(Equation 1) with the following criteria: 


x = log(Srisoow / 93.6) 
y = log(Sr10oow / 93.6) 
x 0.55 


y > 0.05 (3) 


then 13 AGN are selected (out of 22) with a mere five 
contaminants, giving a reliability of 7296. For AGN- 
candidate selection in large surveys, we recommend com- 
bining Equations 1 and 3. However, for small surveys, 
where computational resources are not a concern, we 
recommend template fitting or SED decomposition as 
color selection lacks either completeness or reliability. 


3.5. Number Counts 


We present cumulative number counts of the 10 um 
sources in Figure 8 and Table 4. Overall, there were 
573 galaxies with a >3ø detection at 10 jm (observed 
frame) and 469 10 jum-detected galaxies in the final sam- 
ple that have been fit with templates (these numbers 
differ slightly from Section 2.1 because here we are de- 
scribing sources with a 10 um detection, rather than a 
detection at any wavelength). At the peak redshift dis- 
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Figure 7. MIRI + near-IR color selection of AGN. Symbols 
and lines are the same as in Figure 6. Dust-rich galaxies 
separate from mid-IR weak galaxies, in contrast with Figure 
6. The solid lines show our proposed AGN selection region. 
The main contaminant in this region is dust-rich star-forming 
galaxies at z « 1. 


tribution of the MIRI sources, z — 1.4, the F1000W 
bandpass is probing rest-frame 4 jum emission. 

We compare the 10 um counts with the IRAC 8 um 
counts from Fazio et al. (2004) in the EGS and Boótes 
fields. The EGS counts match the 10jm counts at the 
bright end, although the Boótes counts are higher. In 
the mid-range (10 wJy), the MIRI and IRAC counts also 
disagree. This could possibly be attributed to the much 
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Table 3. Average Spectral Energy Distributions 


Wavelength Ly Unc. Ly Unc. Ly Unc. Ly Unc. 
(um) (W/Hz) (W/Hz) (W/Hz) (W/Hz) (W/Hz) (W/Hz) (W/Hz) (W/Hz) 
log Lig < 9.0L9 log Lig = 9.0 — 9.5Lo log Lig = 9.5 — 10.0L9 | log Lir = 10.0 — 10.5Le 
0.335 | 3.66e20 7.97e19 | 1.24e21 1.57e20 | 2.14e21 2.98e20 | 2.30e21 2.55e20 
0.454 | 6.90e20 1.08e20 | 2.36e21 2.36e20 | 5.96e21 7.42e20 | 6.10e21 6.39e20 
0.544 | 7.91e20 1.21e20 | 2.49e21 5.00e20 | 6.25e21 8.16e20 | 8.19e21 9.12e20 
0.663 | 8.32e20 1.38e20 | 2.98e21 3.06e20 | 6.24e21 5.70e20 | 8.06e21 7.60e20 
0.913 | 1.14e21 1.32e20 | 3.74e21 2.89e20 | 7.78e21 1.06e21 | 9.91e21 1.39e21 
1.330 | 1.20e21 1.41e20 | 4.33e21 3.17e20 | 1.02e22 7.18e20 | 1.44e22 3.56e21 
1.740 | 9.60e20 9.93e19 | 3.69e21 2.16e20 | 9.17e21 7.31e20 | 1.25e22 2.16e21 
2.243 | 8.94e20 4.25e19 | 2.92e21 2.38e20 | 7.14e21 9.80e20 | 1.01e22 1.18e21 
2.750 | 6.63e20 5.02e19 | 1.86e21 1.88e20 | 6.21e21 7.93e20 | 7.99e21 9.35e20 
3.220 | 6.64e20 5.81e19 | 1.93e21 1.82e20 | 4.58e21 4.17e20 | 6.92e21 1.15e21 
3.722 | 6.62e20 6.59e19 | 1.42e21 4.00e20 | 4.67e21 8.32e20 | 6.48e21 1.45e21 
4.462 | 4.67e20 5.32e19 | 1.43e21 9.95e19 | 3.73e21 6.27e20 | 4.65e21 5.52e20 
5.284 | 3.78e20 4.36e19 | 1.48e21 1.29e20 | 4.18e21 6.19e20 | 7.21e21 1.33e21 
5.797 | 4.51e20 6.12019 | 1.63e21 2.54e20 | 5.06e21 5.44e20 | 1.06e22 1.07e21 
6.310 | 4.71e20 5.90e19 | 2.20e21 2.94e20 | 6.89e21 1.01e21 | 1.67e22 1.85e21 
6.784 | 6.38e20 7.38e19 | 3.87e21 4.07e20 | 7.08e21 1.33e21 | 2.05e22 2.08e21 
7.269 | 6.91e20 7.97019 | 4.17e21 3.66e20 | 1.31e22 1.22e21 | 3.00e22 4.30e21 
7.635 | 8.03e20 1.18e20 | 5.03e21 4.15e20 | 1.33e22 2.15e21 | 3.32e22 1.69e21 
8.259 | 8.75e20 1.02e20 | 4.93e21 4.02e20 | 1.24e22 1.08e21 | 2.98e22 1.80e21 
8.683 | 7.23e20 1.21e20 | 3.42e21 6.18e20 | 1.33e22 1.88e21 | 2.68e22 2.66e21 
9.731 | 7.53e20 1.07e20 | 3.71e21 2.88e20 | 8.78e21 6.21020 | 2.11e22 2.85e21 
10.919 | 1.00e21 1.83e20 | 4.37e21 5.27e20 | 1.33e22 1.50e21 | 2.79e22 6.23e21 
12.052 | 1.12e21 1.35e20 | 4.75e21 5.47e20 | 1.39e22 2.67e21 | 3.22e22 3.94e21 
13.908 | 9.27e20 1.91e20 | 2.86e21 2.02e20 | 8.21e21 7.60e20 | 3.80e22 5.43e21 
16.462 | 9.50e20 1.68e20 | 2.76e21 1.85e20 | 1.12e22 1.66e21 | 3.57e22 2.83e21 
log Lin = 10.5 — 11.0Lo log Lir = 11.0 — 11.5Lo log Lin > 11.5Lo 
0.335 | 3.97e21 4.47e20 | 7.31e21 2.75e21 | 5.21e21 2.06e21 
0.454 | 8.55e21 2.11e21 | 8.62e21 1.26e21 | 1.58e22 2.92e21 
0.544 | 1.24e22 3.34e21 | 1.38e22 3.02e21 | 1.94e22 3.28e21 
0.663 | 1.27e22 1.95e21 | 2.08e22 3.08e21 | 3.48e22 6.27e21 
0.913 | 2.00e22 9.34e21 | 2.87e22 6.26e21 | 6.61e22 1.41e22 
1.330 | 1.84e22 5.82e21 | 2.86e22 1.49e21 | 7.97e22 8.70e21 
1.740 | 1.84e22 1.60e21 | 3.24e22 1.63e21 | 9.16622 8.06e21 
2.243 | 1.50e22 1.90e21 | 2.69e22 1.44e21 | 8.44e22 4.86e21 
2.750 | 1.79e22 5.37e21 | 2.27e22 1.59e21 | 6.52e22 9.99e21 
3.220 | 1.15e22 2.24e21 | 2.43e22 1.53e21 | 7.50e22 2.96e21 
3.722 | 1.00e22 1.01e21 | 2.28e22 1.29e21 | 8.07e22 5.04e21 
4.462 | 1.18e22 1.76e21 | 2.07e22 1.32e21 | 8.31e22 6.66e21 
5.284 | 1.37e22 2.26e21 | 3.09e22 4.60e21 | 1.32e23 1.62e22 
5.797 | 2.49e22 4.67e21 | 5.70e22 7.32e21 | 2.53e23 4.62e22 
6.310 | 4.07e22 3.79e21 | 1.00e23 1.34e22 | 4.67e23 2.32e22 
6.784 | 5.54e22 7.34e21 | 1.44e23 2.00e22 | 5.30e23 4.59e22 
7.269 | 8.39e22 9.27e21 | 2.10e23 2.48e22 | 7.79e23 1.06e23 
7.635 | 1.08e23 2.35e22 | 2.52e23 4.26e22 | 1.09e24 7.23e22 
8.259 | 9.49e22 4.55e21 | 2.58e23 3.15e22 | 9.74e23 7.57e22 
8.683 | 9.37e22 1.48e22 | 2.27e23 4.10e22 | 8.25e23 2.20e23 
9.731 | 5.00e22 1.64e22 | 1.50e23 5.35e22 | 5.93e23 1.63e23 
10.919 | 8.87e22 2.47e22 | 1.44e23 2.50e22 | 5.37e23 1.08e23 
12.052 | 7.64e22 1.17e22 | 2.25e23 7.95e22 | 9.58e23 2.87e22 
13.908 | 7.83e22 1.29e22 | 1.92e23 2.91e22 | 6.40e23 9.71e22 
16.462 | 6.51e22 0.00e00 | 1.61e23 9.23e21 | 6.45e23 4.16e22 
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Figure 8. 10 um cumulative number counts for all sources 
with SNR>3 in F1000W (open stars) and the final sample 
(filled stars). Our selection criteria mainly affected sources at 
the faintest F1000W fluxes. We also show the number counts 
for the mid-IR weak galaxies (blue squares), star-forming 
galaxies (f(AGN)mir< 0.5; filled green circles) and AGN 
(f(AGN)urn 0.5; red triangles). AGN predominantly have 
brighter F1000W fluxes. The 8 um number counts from Fazio 
et al. (2004) are shown in pink. The discrepancy between the 
8 um and 104m counts may be due to the large differences 
in field sizes. The number of sources in one MIRI field of 
view is shown on the right axis. 


small area in the MIRI survey compared with the IRAC 
data. A mere 67 sources have Sp1000w = 8.0 — 80.0 uJy, 
which could make the underlying distribution of the 
parent population difficult to distinguish. A large area 
MIRI survey will better determine the true flux distri- 
bution of galaxies. 

We also show the number counts of the mid- 
IR weak, star-forming (f(AGN)mir< 0.5) and AGN 
(f(AGN)mir> 0.5) populations. The faint galaxy popu- 
lation at cosmic noon contains a large fraction of mid-IR 
weak galaxies (Figure 3). Interestingly, both the mid- 
IR weak and dust-rich galaxies have approximately the 
same flux distribution. The mid-IR weak galaxies do 
not dominate the counts at faint fluxes. We find fewer 
than five AGN per MIRI pointing, illustrating that a 
large MIRI survey will be required to build up sufficient 
population statistics of infrared AGN. MIRI is capable 
of finding AGN in lower luminosity galaxies and galaxies 
at higher redshift (Figure 3), but large areas are required 
due to the relative rarity of AGN. 


4. DISCUSSION 
4.1. Where are All the AGN? 


Kirkpatrick et al. (2017a) predicted finding a signifi- 
cant population of faint AGN at cosmic noon by extrap- 
olating down the luminosity function based on the num- 
ber of AGN observed in IR luminous galaxies. The MIRI 
observations have revealed significantly fewer AGN at 
cosmic noon and high redshift than predicted in K17. 
This can be attributed to two reasons: AGN may be 
harder to identify than expected or AGN may by intrin- 
sically low luminosity in lower mass galaxies. Distin- 
guishing between the scenarios will require larger MIRI 
surveys in conjunction with spectroscopic or X-ray sur- 
veys. 

'The most luminous, unobscured AGN exhibit power- 
law emission in the near-IR as well as the mid-IR. 
Heavily obscured AGN, however, may not be visible 
in the near-IR (Hickox & Alexander 2018, and refer- 
ences therein). The emission from the hottest regions of 
the torus will be reabsorbed by subsequent dust layers 
and re-emitted at longer wavelengths. This reprocessing 
weakens the emission at shorter infrared wavelengths. 
'The obscured fraction of AGN rises with increasing red- 
shift (Peca et al. 2023; Yang et al. 2023a), which means 
that relying on rest-frame A < 5 um photometric emis- 
sion will likely miss a significant fraction of AGN beyond 
z > 3, when the MIRI filters no longer cover the mid-IR. 
Heavily obscured AGN may have a visible stellar bump, 
arising from the host galaxy, although power-law emis- 
sion should still be visible in the mid-IR (Lyu & Rieke 
2022). In Figure 7, there are AGN that lie in the mid- 
IR weak region, outside of the proposed AGN selection 
area, indicating the diversity of near-IR emission exhib- 
ited by AGN. AGN with log($1s0ow /53.5) < 0.5 may 
be more obscured. Alternately, these could belong to 
a class of “hot-dust poor" Type 1 AGN, which exhibit 
relatitively weak near-IR emission (Hao et al. 2010; Lyu 
et al. 2017). The cause of the weak torus in these sources 
is unknown, but they have been observed out to z = 6. 

Many galaxies show a mix of power-law AGN emis- 
sion and PAH emission in the mid-IR. In practice, this 
manifests as weaker PAH features (Pope et al. 2008; 
Sajina et al. 2012; Kirkpatrick et al. 2012). Low lu- 
minosity AGN that are accreting at high Eddington 
ratios will have power-law emission from a torus (Ma- 
son et al. 2012), but this emission may be outshone by 
the galaxy’s star forming component. In low luminosity 
galaxies, weak PAH features may be due to swamping 
by an AGN, or they may be intrinsically weak due to 
the galaxy’s dust content, as seen in Figure 5. There- 
fore, at low luminosity, weak PAH features are not a 
reliable means of identifying AGN. Additionally, low lu- 
minosity AGN that accrete at low Eddington ratios may 
lack a torus altogether (Yuan & Narayan 2014). Intrin- 
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Table 4. 10 um Cumulative Number Counts 


S1000w (uy) All Sources Final Sample f(AGN)um&€ 0.5 f(AGN)wir>0.5 mid-IR weak Sources 
0.3 6.9940.30x 108 5.58+0.27x 108 — 2.39 0.18 x 108 Im 2.45 + 0.18 x 108 
0.5 5.6540.27x 108 4.73+0.25x 108 2.08: 0.16 x 108 3.47 0.67 x 107 1.98 + 0.16 x 108 
0.8 4.69+40.25x 108 4.05+0.23x 108 1.85 3-0.15 x 1095. 3.34+0.66 x 107 1.63 + 0.15 x 108 
1.2 3.75+0.22x 108 3.2940.21x 108 1.544 0.14 x 108 2.96+40.62 x 107 1.30 + 0.13 x 108 
2.0 2.88+0.19x 108 2.454018 108 11840.12x108 2.5740.58 x 107 9.00 + 1.08 x 107 
2.9 2.42+0.18x 108 2.034016 108 1.05 20.12 x 108 2.18+40.53 x 107 6.68 + 0.93 x 107 
4.0 1.97+0.16x 108 1.6640.15x108 9.12+1.08x107 1.93+0.50 x 107 4.88 + 0.79 x 107 
5.0 1.68+0.15x 108 1.44: 0.14 x 108 7.842: 1.00 x 107 1.80 0.48 x 107 4.37 + 0.75 x 107 
7.5 1.22:0.13 x 108  1.05:0.11 x 108  5.782:0.86 x 107 1.41 0.43 x 107 3.08 + 0.63 x 107 

15.0 7.3340.97x 107 6.812- 0.94 x 107 — 3.21 0.64 x 107 1.03 0.36 x 107 2.57 + 0.58 x 107 
30.0 4.24+0.74x 10" 3.86 0.70 x 107 1.934-0.50 x 107 5.14 2.57 x 106 1.41 + 0.43 x 107 
50.0 2.96 0.62 x 107  2.702:0.59 x 107 1.16 +0.39 x 107 — 3.86 4 2.23 x 109 1.16 + 0.39 x 107 
75.0 2.06+0.51x 107 1.802: 0.48 x 107 — 7.713: 3.15 x 106 6.43 + 2.87 x 109 
100.0 1.41 +0.43 x 107 1.29 0.41 x 107 5.14 2.57 x 106 3.86 + 2.23 x 106 


sically low luminosity AGN (Lx < 10? erg/s) may be 
difficult or even impossible to identify with photometric 
information alone. 

In both cases (high obscuration or low luminosity), 
there are promising ways forward. Hatcher et al. (2021) 
used mid-IR colors to identify AGN candidates that were 
not X-ray detected in the COSMOS field. The authors 
performed X-ray stacking to obtain the average emis- 
sion of the AGN candidates and found that, indeed, 
as a population, the AGN candidates had X-ray signa- 
tures indicative of hosting a low luminosity AGN. X-ray 
stacking remains a promising way to measure the aver- 
age AGN emission of a population. Specific to JWST, 
nebular line diagnostics will also be useful (Backhaus 
et al. 2022), potentially moreso than looking for broad 
line emission, as low luminosity AGN may lack a broad 
line region (Elitzur et al. 2014). AGN emit multiple 
high ionization nebular lines, including [Nev], [Nevri], 
[Orv] (Cleri et al. 2023a,b; Negus et al. 2023). Mid-IR 
spectroscopy may be the most reliable method for iden- 
tifying obscured or low luminosity AGN (Petric et al. 
2011; Bonato et al. 2017; Stone et al. 2022). 

With photometry alone, headway could possibly be 
made by spatially decomposing the galaxies into an in- 
ner and outer region. For the most massive galaxies, 
this may be possible, but as is clear from Figure 1, 
most MIRI galaxies are small in size. This is likely due 
to the fact that MIRI is probing lower mass galaxies. 
Finally, photometric or spectroscopic identification of 
AGN may be made more reliable through machine learn- 
ing algorithms (Holwerda et al. 2021; Tardugno Poleo 
et al. 2023). Machine learning is capable of combining 
more photometric and color information than is possible 


with simple color selection. Machine learning can also 
take into account redshift information, which is useful 
in identifying the PAH features. 

While identification challenges are a likely culprit for 
not finding more AGN at cosmic noon in our MIRI sam- 
ple, there is also the possibility that bolometrically lu- 
minous AGN do not exist in large quantities in lower 
luminosity (Lig < 10! Lo) galaxies. In isolated, low- 
luminosity galaxies, there may not be a large enough 
gas supply being funnelled to the center of the galaxy to 
fuel black hole growth at high enough Eddington rates 
to be energetically dominant over the star formation and 
therefore detectable. Alternately, but similarly problem- 
atic, is the role that AGN duty cycle may play. AGN 
in low luminosity galaxies may accrete at large enough 
Eddington ratios to be detectable, but they may only do 
so for a short amount of time. If we make the simplify- 
ing assumption that all galaxies host detectable AGN at 
some point in their life, then the left panel of Figure 3 
indicates that the duty cycle of these AGN is < 10% of 
the star-forming time. In this scenario, we will never ob- 
serve large amounts of AGN in low luminosity galaxies 
because they are in their “dormant” phase. 


4.2. PAH Emission 


PAH emission arises at the edges of star-forming re- 
gions, and globally, the integrated PAH emission of a 
galaxy has been shown to correlate with its SFR. The 
link between PAH emission and star formation natu- 
rally leads to a correlation with Lyr, as Lyr is largely 
attributed to dust heating by stars formed within the 
past 100 Myr. PAH emission is frequently parameterized 
by the luminosity of the 6.2 um feature, Le.2. The ra- 
tio Ls.3/ Lig is observed to decrease with increasing Lir. 
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Locally, this decrease is most evident in ULIRGs (Lig > 
107? Lo). The origin of decreasing PAH emission is de- 
bated, but one possibility is that the geometry of star- 
forming regions is changing. In ULIRGs, star forma- 
tion predominantly occurs in a large, obscured nuclear 
starburst. This increases the volume of dust (increas- 
ing Lig) while decreasing the surface area (decreasing 
Lg.2). In contrast, local LIRGs (Lir = 1011 — 10? Lo) 
are more similar to normal star-forming galaxies in that 
they have many star-forming regions spread throughout 
the galaxy. 

Intriguingly, the ratio of L¢.2/Lrp is higher in ULIRGs 
at cosmic noon compared with local ULIRGs, indicat- 
ing changing conditions of star formation (Pope et al. 
2013; Kirkpatrick et al. 2014; McKinney et al. 2020). 
Studies of Lg5/Lig in galaxies beyond the local uni- 
verse are rare, as previously Spitzer/IRS was the only 
instrument capable of measuring Le.2, and only in the 
brightest galaxies. 

We take a first look at estimating L.»/ Lig in galaxies 
at z > 0 using JWST down to Lyr an order of magnitude 
lower than what was possible with Spitzer. The MIRI fil- 
ters cover the 6.2 jum feature in different redshift ranges. 
We identify the redshift ranges for each filter where the 
6.2 um feature is isolated, and the rising continuum be- 
yond 6.3 jum does not dominate the emission. We find 
the following redshift ranges suitable: 


FTTOW: 0.20 < z < 0.28 (4) 
F1000W: 0.58 < z « 0.67 (5) 
F1280W: 1.04 < z < 1.11 (6) 
F1500W: 1.39 < z « 1.47 (7) 
F1800W: 1.87 « z « 1.95 (8) 


We do not include F2100W as we do not have any 
sources in a suitable redshift range. 

For each photometric point, we first estimate how 
much of the broadband photometry is potentially due 
to the 6.2 um emission feature using a local sample with 
mid-IR spectroscopy. The 5MUSES sample was ob- 
served with Spitzer/IRS and have published Ls and 
Lig measurements in Kirkpatrick et al. (2014). That 
work contains 11 star-forming galaxies spanning z = 
0.06 — 0.24, log Lig = 10.79 — 11.63 La, and a range of 
Lo.2/ Lig. The average (standard deviation) of Ls.2/ Lig 
is 0.00456 (0.00181) for the 11 5MUSES galaxies. 

We shift the 11 5MUSES Spitzer/IRS spectra to the 
redshift of each MIRI galaxy and convolve with the ap- 
propriate MIRI transmission filter to calculate v L,, or 
the photometric luminosity in a given MIRI bandpass. 
We use the published Ls» value and calculate L3 /vL, 
for all 11 5MUSES galaxies. We take the mean and 
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Figure 9. Le.2/Līr vs. Lim for the MIRI galaxies falling in 
the redshift ranges in Equation 8. Lo.» has been estimated by 
calculating what percentage of the photometric point could 
be due to the 6.2 um feature using the 5MUSES sample of 
star-forming galaxies from Kirkpatrick et al. (2014, black 
stars). The mid-IR weak galaxies are shown as the open 
symbols, while the filled symbols correspond to the dust rich 
galaxies. Lower luminosity galaxies have higher ratios of 
Le.2/ Lin, in line with local trends. Interestingly, the higher 
redshift sources seems to have higher Le.2/ Lr, which could 
indicate a shift in the star-forming conditions at higher red- 
shift. Spectroscopy is required to accurately measure Le. 
and confirm this trend. 


standard deviation of Lg.s/vL, for the 11 galaxies as 
a scaling factor and associated uncertainty. We then 
multiply the MIRI photometry of each galaxy in our 
sample by this scaling factor to estimate Lg.» [Lo] and 
an uncertainty. Figure 9 shows Lg.2/Lir as a function 
of Lig for the MIRI sources. This is a conservative ap- 
proach that marginalizes over uncertainty in the under- 
lying SED based on what we expect from z « 0.2 galax- 
ies. A caveat to this approach is that Lyr has been de- 
termined by fitting templates to all the MIRI and IRAC 
photometry, which will include the data point used to 
calculate Le.2. While these calculations are somewhat 
independent, spectroscopy is needed to fully uncouple 
them. 

Figure 9 shows the same trend that is observed 
locally-Le.2/Lim decreases with increasing Lig for a 
fixed redshift bin. This decrease starts to happen no- 
ticeably around Lig = 10'° Lo. This is the same lu- 
minosity where the near-IR and mid-IR photometry be- 
gins to be best fit by dusty star-forming templates from 
K15, rather than mid-IR weak templates from CEO1. 
'The conditions of star formation may be significantly 
different in these infrared luminous sources. We show 
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the mid-IR weak sources as the open symbols, and they 
have higher Le.2/Lır ratios. The decrease in Ls.2/ Lin 
with increasing Lig is then driven by a shift from pre- 
dominantly mid-IR weak sources to dust-rich sources. 

Figure 9 also hints at an increase in Lg. 2/ Lig at higher 
redshifts (red triangles compared with the yellow cir- 
cles), similar to what has been measured in sources with 
Spitzer/IRS spectroscopy (Pope et al. 2013). This shift 
could indicate that the geometry and physics of star for- 
mation is changing with increasing lookback time. For 
example, higher Ls »/ Lig ratios at earlier times could be 
a byproduct of the higher gas fractions in high redshift 
galaxies. Indeed, Cortzen et al. (2019) demonstrate that 
PAHs are good tracers of the cold molecular gas from 
z ~ 0— 2, and galaxies at high redshift are increas- 
ingly gas-dominated which would increase Lg. relative 
to the dust-obscured SFR measured by Lig. Alterna- 
tively, McKinney et al. (2020) show that the offset in 
La.2/ Lin for fixed Lig between z ~ 0 and z = 2 galaxies 
disappears when normalizing Lig by the IR size. This 
indicates that the geometry of distributed star formation 
across these galaxies also plays an important role in set- 
ting the Le.2/Lır ratio. Of course, the reader should 
bear in mind two caveats-Lg.» has been estimated from 
photometry and Ljg has been extrapolated from near- 
and mid-IR photometry. Robustly testing how the PAH 
fraction changes with luminosity and redshift will re- 
quire JWST spectroscopic observations and a new far-IR. 
telescope. 


4.3. The Nature of Mid-IR Weak Sources 


We designated sources as mid-IR weak or dust-rich 
based on which set of templates fit them best. Mid-IR 
weak sources have a much larger near-IR/mid-IR emis- 
sion ratio. The mid-IR weak status has a few potential 
explanations. First, these galaxies may be dust-poor, so 
the mid-IR emission overall is weaker than the near-IR 
emission, but the dust grain size distribution (i.e. PAH 
fraction) is similar to the dust-rich galaxies. On the 
other hand, these galaxies may be metal-poor, so that 
the PAH features are weaker. Finally, weaker mid-IR 
emission may be caused by a change in dust heating 
mechanisms. In the local Universe, spatially resolved 
modeling of the mid- and far-IR emission show that the 
temperature of dust in lower mass galaxies is lower due 
to the dusty disk being more extended than the stel- 
lar disk (Trewhella et al. 2000; Dalcanton et al. 2004; 
Xilouris et al. 2004; Holwerda et al. 2009, 2012). The 
weak mid-IR emission is then attributable to most of 
the dust in low mass galaxies being colder. Morpholog- 
ical studies with MIRI and NIRcam may be able probe 


the extent of the dust and stars in z ~ 1 — 2 galaxies 
(Magnelli et al. 2023; Shen et al. 2023). 

The high Ls2/Lig ratios in the mid-IR weak galax- 
ies (Figure 9) hint that decreasing PAH emission (and 
hence low metallicity) is not to blame for the changing 
near-IR/mid-IR ratios. We also test this by fitting the 
average SEDs (Table 3) with the (Draine & Li 2007) 
dust models, in which the spectra are fit with models 
heated by radiation fields of different strengths, for mix- 
tures of amorphous silicate and graphite grains, includ- 
ing varying amounts of PAH particles. We find that all 
but two of the SEDs have the exact same gpar = 4.58 
as measured by these models. The log Lig < 9.0 and 
log Lig = 10.0 — 10.5 SEDs have gpar = 2.50. These 
results imply that the PAH spectrum is quite homoge- 
neous for all of our galaxies, although the PAH emission 
only rises above the 1.6 jum stellar bump for galaxies at 
higher Lig. This is in fact different from how models 
such as CE01 or Rieke et al. (2009) behave, as the PAH 
features in those models become more obscured by a 
dust continuum at higher Lig values. In contrast, the 
K15 templates have much stronger PAH features than 
local templates of the same luminosities. 

If we take Figure 9 at face value, the interpretation of 
the mid-IR weak sources becomes that their overall dust 
emission is weak, but they are not metal poor. The sim- 
ilar gpar values point to similar amounts of very small 
grains regardless of luminosity. The higher Le.2/Lır in 
mid-IR weak galaxies points to a lack of far-IR emis- 
sion, which would rule out the interpretation that their 
dust is colder. Future far-IR observations would pro- 
vide interesting insights into the ISMs of this relatively 
unexplored population. 


4.4. How do MIRI SFRs correlate with Optical SF Rs? 


In dust rich galaxies, infrared star formation rates 
(SFRs) are more reliable than optical SFRs since much 
of the optical emission is attenuated (Kennicutt 1998; 
Kennicutt & Evans 2012; Calzetti et al. 2007; Boquien 
et al. 2016). It is therefore tempting to use IR observa- 
tions to calculate SFRs for all infrared-detected galax- 
ies. However, MIRI observations probe further down 
the luminosity function at cosmic noon than previous 
IR telescopes, and most CEERS MIRI sources are un- 
detected by the available Herschel observations. This 
can introduce significant uncertainty in the calculation 
of Lig as the ratio of far-IR to mid-IR emission has in- 
trinsic scatter and can also depend on factors such as 
metallicity and whether galaxies are on the main se- 
quence (Elbaz et al. 2011). In this subsection, we com- 
pare previously calculated optical SFRs with SFRs es- 
timated by our simplistic template fitting. The purpose 
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Figure 10. Optical vs. IR star formation rates, where sym- 
bol color and shape corresponds to redshift. Mid-IR weak 
galaxies are shown as the open symbols. At lower luminosi- 
ties, the bulk of star formation is unobscured, which is re- 
flected by the optical SFRs being systematically higher than 
the IR SFRs. This is mainly an issue for the mid-IR weak 
galaxies. 


of this section is to provide a first look at how well these 
SFRs compare, and whether different ways of estimat- 
ing SFRs give similar answers. An in-depth analysis of 
MIRI-based SFRs, including new SFR calibrations, will 
be discussed in Roynane et al. (2023, in prep). 

Figure 5 shows how the near-IR and optical emission 
increases at lower infrared luminosities, indicating that 
more of the star formation is unobscured. In the bottom 
panel of Figure 10, we compare IR and optical SFRs for 
our MIRI sample, to examine when the two indicators 
are in disagreement. We calculate the infrared SFRs 
using the equation: 


SFR [Mo /yr] = 1.59 x 101? x Lip [Lo] (9) 


which assumes a Kroupa IMF (Murphy et al. 2011; 
Kroupa 2001). 

'The optical star formation rates come from the CAN- 
DELS catalogs of Stefanon et al. (2017), following meth- 
ods described there (see also Dahlen et al. (2013) and 
Santini et al. (2015)). Briefly, the star formation rates 
are derived by fitting the CANDELS UV/Optical pho- 
tometry in 10 different ways, each fit using a different 
code, priors, grid sampling, and star formation histo- 
ries. The final value is the median from the different fits. 
Optical SFRs are sensitive to the photometric redshift 
of the galaxy. The optical SFRs from Stefanon et al. 
(2017) were computed using photometric redshifts from 
that catalog, whereas for CEERS MIRI galaxies we also 


use photometric redshifts from Finkelstein et al. (2022) 
(see Section 2.2). For this reason, we keep only those 
sources where the photo-z's from Stefanon et al. (2017) 
and Finkelstein et al. (2022) agree within Az « 0.15. 
Similarly, we also remove any sources where the MIRI 
redshift and the optical redshift differ by more than 0.15. 
'The remainder of the sample is shown in Figure 10. 

There is in general a one-to-one correlation between 
the optical and IR SFR, and there is no trend away from 
the one-to-one line with increasing f(AGN)mrr. The 
disagreement is strongest for a handful of mid-IR weak 
galaxies, where the optical SFRs greatly underpredict 
the IR SFRs. At Lig « 10!? L5, the optical SFRs are 
nearly universally higher than IR SFRs. IR SFRs may 
underpredict the true SFR. at low luminosities, when 
the bulk of star formation is unobscured. Interestingly, 
Figure 10 shows that the disagreement between optical 
and IR SFRs depends on whether a galaxy is classified as 
mid-IR weak or dust-rich, rather than on luminosity or 
redshift. For the mid-IR weak galaxies, the optical SFRs 
are higher, while for the dust-rich galaxies, the IR SFRs 
are higher. Dust-rich galaxies at the same redshifts or 
luminosities as mid-IR weak galaxies lie closer to the 
one-to-one line. 

Optical SFRs themselves carry a high degree of un- 
certainty as there are many degeneracies that plague 
this portion of a galaxy's SED, including stellar ages 
and dust attenuation (Leja et al. 2019; Bell 2003). This 
comparison between SFR indicators is for optical- and 
IR-based SFRs calculated with template-fitting to a few 
photometric data points. Roynane et al. 2023 (in prep.) 
compares mid-IR-based SFRs with UV-based SFRs and 
finds a much stronger agreement. The most accurate es- 
timations of SFRs will likely come from combining either 
UV or spectroscopic indicators, such as Ha, with mid- 
IR estimators (Murphy et al. 2011; Kennicutt & Evans 
2012). 


5. CONCLUSIONS 


We have examined the demographics of the CEERS 
MIRI survey, which comprises four pointings in the EGS 
field covering the F770W-F2100W filters. We summa- 
rize our findings below. 


e Comparison with the Spitzer 24 um popula- 
tion: MIRI sources peak at z 1.3 and lack many 
high-redshift detections seen in the MIPS 24 um 
population. This is due both to MIRI’s smaller 
field of view and to MIRI’s lower wavelength filters 
(which are more sensitive), no longer tracing dust 
beyond z > 3. A large area MIRI survey will be re- 
quired to build up a large high-redshift population. 
However, in much shorter integration times, MIRI 
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probes an order of magnitude fainter than MIPS 
due to the increased sensitivity of JWST. MIRI is 
best at identifying galaxies with Lig < 101? Lo at 
cosmic noon (z = 1 — 2). 


e Mid-IR weak galaxies: Until now, with the im- 
proved mid-IR coverage and sensitivity of JWST, 
it was impossible to know exactly what the mid- 
IR emission of Lig < 101° Lo galaxies at cosmic 
noon looked like. We fit the MIRI and IRAC 
4.5 um photometry with a suite of dusty star form- 
ing templates, AGN templates, and templates cre- 
ated from local low-luminosity galaxies. We find 
that inclusion of the 4.5 um data point is essential 
to discriminate between AGN and “mid-IR weak” 
galaxies, as their mid-IR emission can look very 
similar. It is only by comparing the mid-IR emis- 
sion with rest-frame near-IR emission that one can 
distinguish between a mid-IR power-law due to an 
AGN torus and weak mid-IR emission due to an 
intrinsically low-luminosity galaxy. We find that 
mid-IR weak galaxies dominate the MIRI popula- 
tion at Lig < 10'° Lo. In these galaxies, the SED 
is dominated by near-IR emission, indicating that 
the star-formation is predominantly unobscured. 
For that reason, we recommend that infrared SFRs 
should be used with extreme caution. 


e AGN identification with colors: The AGN 
fraction increases with Lig, which may reflect 


an increasing accretion efficiency, or it may re- 
flect the challenges in separating AGN from mid- 
IR weak galaxies at lower luminosities. Spec- 
troscopic surveys are necessary to distinguish be- 
tween the two scenarios. We find that AGN can 
be reliably identified by combining two color se- 
lections: Sigoow/Si000ow V- S1280w /S1000w and 
Sigoow/S3.6 V. Siooow/53.6- The 10 um num- 
ber counts reveal that AGN are extremely rare (a 
handful appear in each MIRI pointing), so large 
area MIRI surveys are required to build up suffi- 
cient samples of mid-IR AGN to study supermas- 
sive black hole growth at cosmic noon. 


MIRI is revealing many interesting properties of galax- 
ies at cosmic noon, but JWST results are hampered by 
it’s small field of view and it’s lack of longer wavelength 
filters. All the results in this paper will be strengthened 
with the upcoming MEGA survey (PI Kirkpatrick, Cy- 
cle 2 GO-3794) which will observe the EGS field with an 
additional 26 MIRI pointings. A new far-IR telescope 
that matches JWS'T's sensitivity would also allow for 
the unambiguous measurement of Lir. 

We thank the anonymous referee for the time and 
energy they spent improving the quality of this pa- 
per. We gratefully acknowledge support from NASA 
grants JWST-ERS-01345 and JWST-AR-02446. This 
work was based on observations made with the 
NASA/ESA/CSA James Webb Space Telescope. 
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